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T^ , Abstract 

OO ■ 

T^lj- ■ We review the resummation of threshold logarithms for heavy 

CN , quark, dijet, direct photon, and W boson production cross sections 

^sT I in hadronic collisions. Beyond leading logarithms the resummed cross 

Q^ ■ section is sensitive to the color exchange in the hard scattering. The 

(— I , resummation is formulated at next-to-leading logarithmic or higher ac- 

0^1 curacy in terms of anomalous dimension matrices which describe the 

o ' factorization of soft gluons from the hard scattering. We give results 

1^ . for the soft anomalous dimension matrices at one loop for the full range 

of partonic subprocesses involved in heavy quark, dijet, direct photon, 

and W boson production. We discuss the general diagonalization pro- 

^ ■ cedure that can be implemented for the calculation of the resummed 

cross sections, and we give numerical results for top quark production 

at the Fermilab Tevatron. We also present analytical results for the 

one- and two-loop expansions of the resummed cross sections. 
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1 Introduction 



The calculation of a large number of hadronic cross sections in perturbative 
Quantum Chromodynamics (QCD) relies on factorization theorems |jl|, H, ^, ^ 
which separate the short- distance hard scattering, calculable in perturbation 
theory, from universal parton distribution functions that are determined from 
experiment. Factorization introduces a scale /x which separates the short- and 
long-distance physics; the factorization scale is often taken to be the same 
as the renormalization scale, though in principle these two quantities are in- 
dependent. Leading order calculations for the cross sections exhibit a strong 
dependence on /i, but next-to-leading order (NLO) calculations diminish this 
scale dependence, and sensitivity to even higher orders is normally exhibited 
by the variation of the scale. NLO results exist for a wide range of processes 
including Drell-Yan p, direct-photon |Q, electroweak boson [|^, Higgs [§], 
heavy quark f^, and jet production |l^. Results exist at two-loops or higher 
for some processes such as Drell-Yan production |11[], e~^e 
and deep inelastic scattering 



annihilation 12 



13 



Calculations to even higher orders are 
formidable and in general one cannot make more precise predictions. Near 
threshold for the production of the observed final state, however, there are 
large logarithmic corrections which arise from soft gluon emission and which 
can be resummed to all orders in perturbation theory |T^. At threshold there 
is no phase space left for the emission of gluonic radiation and the cancella- 
tion of infrared divergences between virtual and real diagrams is incomplete. 
The resultant singular distributions were first resummed for the Drell-Yan 
process [13, W\] this resummation is related to an earlier study of the trans- 
verse momentum distribution of hadron pairs in e+e~ annihilation Jl^. Here 
we will review the threshold resummation of Sudakov logarithms in heavy 
quark, dijet, direct photon, and W boson production in hadronic collisions, 
giving special attention to the resummation of next-to-leading logarithms. 
The resummation of the leading theshold logarithms in heavy quark pro- 



duction cross sections 
|5 



m 



m 



distributions 



2^ and inclusive differential 
I9| pl[] is based on the analogy with the Drell-Yan process 
since the leading logarithms (LL) are universal. The top quark has finally 
been found at the Fermilab Tevatron ^, ^] (for a review see p8| , p9|). The 
top's mass and production cross section are of great interest, and the pre- 
cision of their measurements will increase in future runs, so it is important 
to have good theoretical predictions for the cross section in terms of the top 



mass; threshold effects are expected to play a significant role in these pre- 
dictions. Threshold resummation will also be of relevance to bottom quark 
production at the HERA-B experiment at DESY. 

A study of subleading logarithms in heavy quark production near thresh- 



old and an attempt to resum them was first presented in Ref. []19|. The 
resummation of next-to-leading logarithms (NLL) for QCD hard scattering 
and heavy quark production, in particular, was first put on a solid theoret- 
ical framework in |21|, ^, |31|. It was found that at NLL one has to take 
into account the color exchange in the hard scattering. The resummation 
of the cross section is achieved by refactorizing the cross section into hard 
components associated with the hard scattering, center-of-mass distributions 



for the colliding partons, and a soft function that describes coherent P2 



non- 



collinear soft gluon emission. The soft function is a matrix in color space and 
it satisfies a renormalization group equation in terms of soft anomalous di- 
mensions. The resummation is then given in terms of the angular- dependent 
eigenvectors and eigenvalues of the anomalous dimension matrix. Applica- 
tions to top and bottom quark production at a fixed center-of-mass scattering 
angle, 9 = 90° (where the anomalous dimension matrix is diagonal), were dis- 
cussed in Ref. |Q. More recently the full angle-integrated cross section was 
presented and numerical results given for top quark production at the Teva- 
tron [0]. A different formalism [^] which avoids the angular dependence in 
the resummed corrections has appeared recently for the heavy quark total 
cross section only (and also for direct photon production [Q). The formal- 
ism of Refs. 1^, ^ ^ and its extensions that we review in this paper is 
more general since it keeps an explicit angular dependence; it can thus also 
be applied easily to differential distributions p7| . 

Natural extensions of the NLL resummation formalism to the hadropro- 
duction of dijets have been presented in Refs. [^, ^. In addition, rapidity 
gaps in high-pT dijet events have been studied by resumming logarithms 
of the energy flow into the central rapidity interval [^|. The resumma- 



tion formalism has also been applied to single-particle inclusive cross sec- 
tions 1^ (including direct photon production) and the electroproduction of 
heavy quarks ||4^. Extensions to W + jet production are also straightfor- 
ward Q. 

The NLL resummation formalism that we review here concerns hadronic 



processes of the form 

ha{Pa) + h{pb) ^ T{Q\ y,x) + X, (1.1) 

where ha-, hb, are incoming hadrons and T{Q'^,y,x) denotes the hadronic 
final state, such as a pair of heavy quarks or high-p^ jets, produced with 
large invariant mass Q and rapidity y. The variable x represents the internal 
structure of the final state, for example the angle 6 between the direction of 
the produced heavy quark, or jet, and the beam axis, or the rapidity interval 
between the two final-state jets. The cross section for such processes may be 
written in factorized form as a convolution of the perturbatively calculable 
hard-scattering partonic cross section, i^/„/^, with parton distributions (pf^/h, 
at factorization scale /i, for parton /j carrying a momentum fraction Xi of 
hadron h. Thus, we have: 

dah,h^^T{S,Q^,y,x) s^ r dxa dxb . , 2n , / 2^ 

x^/«/.frv^'^'X,^,«s(/i')) , (1.2) 

where S is the center-of-mass (cm.) energy squared of the incoming hadrons 
and we sum over all relevant partons /a, /{,. Although the physical cross 
section cannot depend on the choice of factorization scheme (usually DIS or 
MS) and factorization scale /x, at any fixed order there is such a dependence. 
Hj^j^^ includes singular distributions which arise, as we mentioned before, 
from incomplete cancellations between cross sections with gluon emission 
and with virtual gluon corrections. These are "plus" distributions, of the 
form [ln™(l — -2)/(l — -2)] + , m < 2n — 1 at nth order in the perturbative 
expansion, which are singular for z = 1, where 

- ^ ^, (1.3) 



with s = XaXbS the invariant mass squared of the partons that initiate the 
hard scattering. We shall refer to 2; = 1 as "partonic threshold" or the 
"elastic limit". We note that by partonic threshold, we mean that the cm. 
total energy of the incoming partons is just enough to produce a fixed final 
state, such as a pair of heavy quarks or jets; thus, heavy quarks, for example, 
are not necessarily produced at rest. 



The exponentiation of singular distributions is derived in terms of mo- 
ments of the cross section with respect to 

r=^. (1.4) 

S ^ ' 



Under moments, the convolution in Eq. ( |1.2| ) becomes a product of moments 
of the parton distributions and the partonic cross section. The exponenti- 
ated cross section in moment space must be inverted back to momentum 
space to derive the physical cross section. There are a number of differ- 
ent prescriptions (and an ensuing debate in the literature) for implementing 
this inversion and avoiding the Landau pole that have been proposed for the 
Drell-Yan process P3, H3|, M and heavy quark production [f^, E^, O. These 



prescriptions do not of course exhaust all possibilities. Here, we will mostly 
be concerned with the resummed cross section in moment space. 

In Section 2 we present the resummation formalism for heavy quark pro- 
duction. In Sections 3 and 4 we give explicit results for the soft anomalous 
dimension matrices in the channels qq — »■ QQ and gg —>■ QQ, relevant to 
heavy quark production. Some technical details of the calculation are given 
in the Appendix. We also give one- and two-loop expansions of the resummed 
cross section. The anomalous dimension matrices are in general not diago- 
nal. In Section 5 we discuss the diagonalization procedure and present some 
numerical results for top quark production at the Fermilab Tevatron. In Sec- 
tion 6 we discuss threshold resummation for dijet production and consider 
the additional complications due to the final state jets. In Sections 7, 8, 
and 9, we give results for the anomalous dimension matrices for the many 
subprocesses relevant to dijet production. In Section 10 we briefly discuss 
resummation for direct photon and W boson production in the context of 
single-particle inclusive kinematics. We conclude with a summary. 

2 Threshold resummation for heavy quark 
production 

In this section we review the general formalism for the resummation of thresh- 
old singularities for heavy quark production. We first write the factorized 
form of the cross section and identify singular distributions in it near thresh- 
old. Then we refactorize the cross section into functions associated with 



gluons collinear to the incoming quarks, non-collinear soft gluons, and the 
hard scattering. Resumniation follows from the renormalization properties 
of these functions and is given in terms of anomalous dimension matrices for 
each partonic subprocess involved. 

2.1 Factorized cross section 

We consider the production of a pair of heavy quarks of momenta pi, p2, in 
collisions of incoming hadrons ha and h^ with momenta Pa and pb, 

haiPa) + h{pb) -^ Q{Pl) + Q{P2) + X , (2.1) 

with total rapidity y and scattering angle 6 in the pair center-of-mass frame. 
The production cross section can be written in a factorized form as 

d(^hah,^QQ{S,Q\y,0) ^ f dXa dXb , 2^, , 2n 

ff=gg,99 

xHjf(-^,y,e,^,as{fi')] , (2.2) 

where S = {pa + PbY, Q^ = (pi + ^2)^, and we sum over the two main 
production partonic subprocesses, qq -^ QQ and gg — > QQ. The short- 
distance hard scattering H^j is a smooth function only away from the edges 
of partonic phase space. The parton distribution functions are given in terms 
of the partonic momentum fractions and the factorization scale fi, which 
separates long-distance from short-distance physics. 

By using the observation of [^ that we can treat the total rapidity of 
the heavy quark pair as a constant, equal to its value at threshold, we can 
rewrite the above cross section as 



d(Th^h,-.QQiS,Q^y,e) _ ^ /-l fdXadXb 2n 

dQ^dyd COS 9 - „-^ Jr'^'J^^'^'^'^^''^'^' 

^ ^ ff=qq,99 



Q 



^^ff^QQ\l-z,^,e,as{^^)\ , (2.3) 



where in addition we have introduced an exphcit integration over z = Q'^/s, 
and a simphfied hard scattering function a. 

In general, a includes "plus" distributions with respect to 1 — z, with 
singularities at nth order in a^ of the type 
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in"*ri 



m < 2n — 1 . 



(2.4) 



These "plus" distributions are defined by their integrals with any smooth 
functions J-'iz) (such as parton distribution functions) by 



dz 
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(2.5) 



All distributions of this kind have been resummed for the Drell-Yan produc- 
tion cross section at leading and nonleading logarithms |T^, |1^] . Here we will 
review the more recent work on resummation for heavy quark production 
including next-to leading logarithms pi], BO, RT . 



In order to calculate the hard-scattering function a ^j^qq^ we consider the 
infrared regularized parton-parton scattering cross section, which factorizes 
as the hadronic cross section. After we integrate over rapidity, we have 



d(Tff_QQ{S,Q'^,9,e) 



dQ"^ dcosO 



X S \ z 



dz 



uXa ClXfy , 2 \ 1 r 2 \ 



a 



ff- 



l-z,^,^,a,(/i2)l . (2.6) 



The argument e represents the universal coUinear singularities. We note that 
the leading power as z ^ 1 comes entirely from the flavor diagonal parton 
distributions 0///(a;a,/i^, e) and (pju^xi^, jj?,e). 

If we take Mellin transforms of the above equation, the convolution be- 
comes a simple product of the moments of the parton distributions and the 
hard scattering function a: 

dQ"^ d cos 6 

-- 4>f/f{N, fi\ e) 4>f/f{N, /i', e) %-.qq(A^, Q/i^, 6, a^ifi^)) , (2.7) 



with the moments given by cr(A^) = /q dz z^^^a{z), (f){N) = /q dx x^~^(f){x). 
We then factorize the initial-state colhnear divergences into the parton distri- 
bution functions, expanded to the same order in a^ as the partonic cross sec- 
tion, and thus obtain the perturbative expansion for the infrared-safe hard- 
scattering function, a. 

We note that under moments divergent distributions in 1 — z produce 
powers of In N: 

/' dz z""-' \^]n}:zA\ = (zirii i^m+1 j^^^ /j^^-i j^\ ^2.8) 
Jo [ 1- z \^ m + 1 V J ^^ 

The hard scattering function a still has sensitivity to soft-gluon dynamics 
through its 1 — 2; dependence (or the A^ dependence of its moments). We 
may now refactorize (moments of) the cross section into A^-independent hard 
components Hjl, which describe the truly short-distance hard-scattering, 
center-of-mass distributions tp, associated with gluons colhnear to the in- 
coming partons, and a soft gluon function 5*^/ associated with non-collinear 
soft gluons. / and L are color indices that describe the color structure of the 
hard scattering. Then we write the refactorized cross section as pi 

dQ"^ d cos u iL \l^ J 



^,^,«.(/.^))^//;(iV,^,eU,-/,-(iV,^, 



X ~Su [ ^, e, a,(/x2) ) ^^/^ I iV, ^, e 1 i^j/j ( iV, -, e ) + 0{l/N) . (2.9) 



This factorization is illustrated in Fig. 1, for the process 

f{Pa) + Rpb) - Q{Pl) + Q{P2) , (2.10) 

in which // represents a pair of light quarks or gluons that annihilate or 
fuse, respectively, into a heavy quark pair. 

The hard-scattering function takes contributions from the amplitude and 
its complex conjugate, 

HjL p, e, a,{^x')\ = h\ p, e, «,(/i^)^ hi f ^, e, as{ii')\ . (2.11) 

The center-of-mass distribution functions ip absorb the universal colhnear 
singularities associated with the initial-state partons in the refactorized cross 




»LI 



(a) 







"^^m 




(b) 

Fig. 1. (a) Factorization for heavy quark production near partonic threshold, 
(b) The soft-gluon function Su, in which the vertices c/, c^ hnk ordered 
exponentials. 
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section. They differ from standard light-cone parton distributions by being 
defined at fixed energy, rather than hght-hke momentum fraction. We define 
them by analogy to the light-cone parton distributions via the matrix 
elements, evaluated in n ■ A = axial gauge in the partonic cm. frame, 

1 r°° ^ 1 

V',/,(x,2po//i,e) = TT^^ dyoe-'-P<^y"{qipMyo,0)-v-^qiOMp)), 



27r23/2y_oo ^ N-.v/'/i^v^u, ^2 



ipg/gix, 2po//i, e) 



ipg/g{x,2po/fx,e) 



27r23/2 

X r t/|/oe-"'P°^»(g(p)|Tr^"^ 



1 



-v--fq{yo,0)q{Q) 



kip)) , 



27r23/2p+ 

/oo 
dyo e-'^P'y^'{g{p)\v,F'^^{yo,0)v,F''^iOMp)) , 

(2.12) 

where the vector v is light-like in the opposite direction from p^, and the 
argument e denotes the coUinear singularities that ip absorbs. The moments 
of ip can then be written as products of moments of the light-cone parton 
distributions and an infrared safe function. 

2.2 Eikonal cross section 

The soft function 5"^/ represents the coupling of soft gluons to the partons 
in the hard scattering process. This coupling may be described by ordered 
exponentials [^, or eikonal or Wilson lines, written as |]38| , |39 



^^/\\2,\i;x)=Pexpl-ig f ' d\p-A^f\XP + x)\ , (2.13) 



where the gauge field A^'^' is a matrix in the representation of the parton 
flavor / of the gauge group SU{3), [3 is the four- velocity of the corresponding 
parton, and P is an operator that orders group products in the same sense 
as ordering in the variable A. 

The color tensors of the hard scattering connect together the eikonal lines 
to which soft gluons couple. We can construct an eikonal operator describing 



11 



soft-gluon emission as 



<l>i{'Hoo,0:x 



.if- 



'P2 



, U, X)c2,d2 ^P 



<l>^-^;^(oo,0;x),,,,, 



X 



■Ma) 



Ub), 



ci),,d.,dA *t (0' -^5 ^)<i.,c.*K (0' -°o; ^)^. 



'b,Cb ' 



(2.14) 

with c/ a color tensor. Then, we may write a dimensionless eikonal cross 
section, which describes the emission of soft gluons by the eikonal lines as 



aif)(^^,^,a.(/x^),e) = E'^(^-^(0) 



X (0|f [K(0))|,^^] 10 m [wj{0)^ |0) , (2.15) 

where ^ is a set of intermediate states which contribute to the weight, w. 
The moments of o"'''^*^ can be factorized, like the moments of the full cross 
section, into moments of the soft gluon function, 5", times moments of jet 
functions, jm, analogous to the ip^s, which absorb the collinear divergences 
of the incoming eikonal lines. Hence S is free of these collinear divergences. 
Thus, we may write 



(eik) 



'Q_ 



,9,as{fi^),e 



S 



LI 






,9,as{fi^ 



xjil^' 



Q 



,as{fi^ 



~Afb) 

JIN 



Q 



,«,(/i2),e , (2.16) 



where the incoming eikonal jet functions are given by products of eikonal 
lines as [^, ^ 



Jin 



{/,) (wiQ 



/^ 



a,{fi^),e 



27r 



dyo e 



-iWiQyo 



X (0| Trj f[^^/f\0, -oo; y)] T[^^/:\o, -oo; 0)] } |0) , (2.17) 
with y'-' = (yo, 0) a vector at the spatial origin. 



2.3 Resummation 

Now, comparing Eqs. ( |2.7| ) and ( |2.9|) , we see that the moments of the heavy- 
quark production cross section are given by 

,n 2 



^ff^QQ 



(AT) 






^HjL —,9,as{iJ,' 

IL \f^ 
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X S 



LI 



'Q_ 



,9,as{fi^ 



(2.18) 



where // denotes qq or gg, the sum over / and L is over color tensors, and 
we have used the fact that ipg/g = ipq/q and (pg/g = (pq/q- Each of the factors 
in Eq. (|2.18|) is gauge and factorization scale dependent. The constraint 
that the product of these factors is independent of the choice of gauge and 
factorization scale results in the exponentiation of logarithms of A^ in ip/ ip 



and Sli [19, 50 1. We now proceed to discuss exponentiation for each factor. 



2.3.1 Parton distributions 



The first factor, (^////0///)^, in Eq. (|2.18|) is "universal" between elec- 
troweak and QCD-induced hard processes, and was computed first with f = q 
for the Drell-Yan cross section |TH]. First, let's present the resummed expres- 



(2.19) 



sion for the ratio ip/(j), with ^ = Q. We have 






%)(a,(Q2)) exp[i?(^)(iV,Q2 



where 



E^f^ (n, Q 



,N-1 



dz- 



\-z 
+ 5(/)k((l-^)™^Q2^ 



'{\-zf A L 



+ i.« 



a.((l-2)2Q2) (2.20) 



and R(^f){as) is an iV-independent function of the coupling, which can be nor- 
malized to unity at zeroth order. The parameter ms and the resummed co- 
efficients B^-f^ depend on the factorization scheme. This scheme-dependence 
must be compensated for by differences in the parton distributions them- 
selves. In the DIS and MS factorization schemes we have ms = 1 and 
ms = 0, respectively. 

^(/)^^(/) and u^f^ are finite functions of their arguments and below we 
give expressions for them needed at next-to-leading order accuracy in InA^. 
We have ^ 



^(/)( 



a. 




(2.21) 



13 



with Cf = Cp = (A^c "~ l)/(2^c) for an incoming quark, and Cf = Ca = Nc 
for an incoming gluon, with Nc the number of colors. Also 



where n/ is the number of quark flavors [^. B^^^ is given for quarks in the 
DIS scheme by 

B(^)(as) = -^Cp^, (2.23) 

4 vr 

while it vanishes in the MS scheme for quarks and gluons. Note that the 
DIS scheme is normally applied to quarks only, but extended definitions 
for gluons are possible [^. Finally, the lowest-order approximation to the 



scheme-independent u^-'' is |^ 



a 

IT 



(/) = 2Cf—. (2.24) 



The above results were ior fi = Q. To change the scale fi, we need the 
renormalization group behavior of the parton distributions ip and 0. 

The distribution ip, and each of its moments, renormalizes multiplica- 
tively, because it is the matrix element of a product of renormalized opera- 
tors. Thus, it obeys the renormalization group equation 

^^hAl^^RlM = 27,(«.(/i^)) ^f/AN, Q/;., e) , (2.25) 

where 7j is the anomalous dimension of the field of flavor /, which is inde- 
pendent of A^. 

The evolution of the light-cone distribution {/)/// with the factorization 
scale fi depends on the factorization scheme that we choose, such as MS or 
DIS. The simplest case is the MS factorization scheme. In this scheme, the 
moments of obey the renormalization group equation 



e 



^d^f^AN^ = 27,KA^,«.(/z^)) <Pf/AN,f,'\e) , (2.26) 

where 7// is the anomalous dimension of the color-diagonal splitting function 
for flavor / |^ . Since only color-diagonal splitting functions are singular as 



X — »• 1, only the flavor-diagonal evolution survives in the large A^ limit. 
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Then the generahzation of Eq. ( p. 191 ) is 
Vi///(iV,Q//i,e) 



0///(^>/^^' 



%)(a,(/i2)) exp[E^f\N,Q') 



Q dfi' 



X exp -2 / -^ 7/(«.(/^'')) - lff{N, as{^")) 
I J fi /i '- 



(2.27) 



2.3.2 Renormalization of the hard and soft functions 

Next, we discuss resummation for the soft function. The soft matrix 5*^/ 
depends on A^ through the ratio Q/{Nn), and it requires renormahzation 
as a composite operator. Its A^-dependence, then, can be resummed by 
renormahzation group analysis |Q ^ Q 0. However, the product Hj^Slj 
of the soft function and the hard factors needs no overall renormalization, 
because the UV divergences of Sli are balanced by construction by those of 
HjL. Thus, we have pl[ 



H 



(0) 
IL 



s 



(0) 



i=a,b 

4- 



-1 



Zs' 



IC 



H, 



CD 



DL 



'LI — {Zs)lbSbaZs,ai-, (2.28) 

where H'^^^ and S^^^ denote the unrenormalized quantities, Zi is the renormal- 
ization constant of the iih incoming partonic field external to hi, and Zs^li 
is a matrix of renormalization constants, which describe the renormalization 
of the soft function, including mixing of color structures. Zs is defined to 
include the wave function renormalization necessary for the outgoing eikonal 
lines that represent the heavy quarks. 

From Eq. ( 2.28|) , we see that the soft function Sli satisfies the renormal- 



ization group equation 31 



!4^^<«^£ 



Sli — —(Xs)lbSbi — Sla(Xs)ai , 



(2.29) 



where Ts is an anomalous dimension matrix that is calculated by explicit 
renormalization of the soft function. In a minimal subtraction renormal- 
ization scheme and with e = 4 — n, where n is the number of space-time 
dimensions, the matrix of anomalous dimensions at one loop is given by 

9 d 



^sig) 



2dg 



Res.^oZsig, e) 



(2.30) 
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Explicit results for the soft anomalous dimension matrices Ts for the partonic 
subprocesses involved in heavy quark production will be presented in the next 
two sections. 

The renormalization group equation ( p.29| ) is, in general, a matrix equa- 
tion. Its solution is 

'Q 



TTiH\-,e,as{n^)\ Pexp 



(2.31) 



xS{l,9,ajQyN')) Pexp[/''^''^r5(a.(^'^))l| 

where the trace is taken in color space, so that 

Tt[HS] = HjlSli = hjhlSu . (2.32) 

At lowest order, Su = Tr [c^^c/]. The symbols P and P refer to path ordering 
in the same sense as the integration variable /i' and against it, respectively. 



2.3.3 Resummed cross section 



Using Eqs. ( p.l8| ), ( p.27|) , and ( p.31|) , the resummed heavy quark cross section 
in moment space is then 

rQ dfi' 



-2 



'Q 



^Ji jj' 



7/,(as(/i' ))-7/,/,(iV,a.(/i' )) 



xTi\H{ —,e,asiiJ.^)\ Pexp 



/ ^rUa,(/i'2) 



S (l, 9, asiQ'/N')) Pexp \ f^'' % Ts («.(/i'^))] | 



(2.33) 

At the level of leading logarithms of A^ in the soft gluon function S, and 
therefore at next-to-leading logarithm of A^ in the cross section as a whole, 
we may simplify this result by choosing a color basis in which the anomalous 
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dimension matrix Ts is diagonal, with eigenvalues A/ for each basis color 
tensor labelled by /. Then, we have 



S 



LI 



'Q_ 



9,as{n^) 



X exp 



SLlU,0,aslj;^ 



Q/N fi 



. (2.34) 



Thus, in a diagonal basis, and with fi = Q and i?(/) normalized to unity, 
we can rewrite the resummed cross section in a simplified form as 



a 



ff^QQ 



(AT) 



Q 



HrL I ^, e, as{l2^) j Sli (l, 9, asiQ'/N^) 
'EW\N,e,Q' 



X exp 



(2.35) 



where the exponent is 



1 ^Jv-i _ 1 



i^-^rs d\ 



iff) 



9Y"[as{XQ' 



1 — z [J{i-z)2 A 
+ gPMil - zr^Q')] + gi'^^'^'^asiil - zrQ'),9]} , 

(2.36) 



with the functions gi, g2, and gs defined by 



9i 



93 



(ILJf) 



-\j - A* + iz/(^) + -u^f^ 



(2.37) 



In the next two sections we present the soft anomalous dimension matri- 
ces for heavy quark production through light quark annihilation and gluon 
fusion, and we give one- and two-loop expansions of the resummed cross 
section. 



3 Soft anomalous dimension matrix for 
the process qq — > QQ 
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3.1 Eikonal diagrams and Ts for qq -^ QQ 

We begin with the soft anomalous dimension matrix for heavy quark pro- 
duction through hght quark annihilation, 

q{Pa, Ta) + q{pb, n) -^ Q{pi, ri) + Q{p2, ra) , (3.1) 

where the p^'s and r^'s denote momenta and colors of the partons in the 
process. The anomalous dimension matrix for this process is 2 x 2, since the 
color indices I and L range over two values. We introduce the Mandelstam 
invariants 

s = iPa+Pbf, h = {pa - Pif - m^ , ui = {pb - Pif - ni^ , (3.2) 

with m the heavy quark mass, which satisfy s + ti + mi = at partonic 
threshold. The UV divergent contribution to S^j is the sum of graphs in 
Figs. 2(a) and 2(b). We give details of the calculation in the Appendix. 

In our calculations we use Feynman rules for eikonal diagrams in axial 
gauge as shown in Fig. 3 (resummation can be performed in a covariant gauge 
as well |]S^). We define dimensionless vectors vf by pf = Qff/-\/2, which 
obey vf = for the light incoming quarks and vf = 2m? /Q'^ for the outgoing 
heavy quarks. The propagator for a quark, antiquark or gluon eikonal line is 
then given by 

. \. , (3.3) 

where 5 = +1(— 1) for the momentum q flowing in the same (opposite) di- 
rection as the dimensionless vector v. The interaction vertex for a quark or 
antiquark eikonal line is 

-WsiT^F)baV^A, (3.4) 

with A = +1(— 1) for a quark (antiquark). The Tp are the generators of 
SU{3) in the fundamental representation. 

For the determination of Fg an appropriate choice of color basis has to 
be made, e.g. singlet exchange in the s- and w-channels, or singlet and octet 
exchange in the s-channel. Here we give the result in a general axial gauge 
in the s-channel singlet-octet basis: 

1 1 

^2 = Coctet = K-'- F/rbraK-'- F)r2ri ~ ~ ^TT^ Csinglet + :^'Jrar2^ri,ri • (3.5) 
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iT(riT(riT(rff'(r(rf(r6'iy6'(r6'^6'(r6'(r6'(r6's^ 





^6'iym 



m^^^^^' 



(a) 





(b) 

Fig. 2. One-loop corrections to Su for partonic subprocesses in heavy quark 
or dijet production: (a) eikonal vertex corrections; (b) eikonal self-energy 
graphs for heavy quark production. 



v7(-v-k+i8) 



-v^/(-v • k+i 8) 
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-v7(v-k+i8) 




v^/(vk+i8) 



Fig. 3. Feynman rules for eikonal lines representing quarks, in the process 
qq -^ QQ. The gluon momentum flows out of the eikonal lines. Group 
matrices at the vertices are the same as for quark lines. 



The results of our calculation are [0, 
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with 



I'll ri2 



Til 

ri2 



-Cf [Lf3 + ln(20v^) + ni] 

IT Ca ^ V tl 

a. 



7r 



Cf 



\tiJ 


- \n{2^yualyb) -Lp-Tii 


31nf"0 

Vti / 


— In + L/3 + TT? 

\tiUi) 



-^ 



where Lp is the velocity-dependent eikonal function 

I -2m' Is /, 1-/3 



^/3 



2mV. / ^ 1 



/3 



1 + /5 



+ TTi 



(3.6) 



(3.7) 



(3.8) 
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with [3 = J I — Am? / s. Note that here we haven't absorbed the function z/^-^-*, 
Eq. ( p.24|) , in the results for Fg, as was done in Refs. [^, ^ 0. The gauge 
dependence of the incoming eikonal hues is given in terms of 



v.; = 



{Vi ■ nf 



\n\ 



(3.9) 



This gauge dependence cancels against corresponding terms in the parton 
distributions. The gauge dependence of the outgoing heavy quarks is can- 
celled by the inclusion of the self-energy diagrams in Fig. 2(b) (see also the 
discussion in the Appendix). 

Vs is diagonalized in this singlet-octet basis at absolute threshold, /? = 0, 
and also for arbitrary 13 when the parton-parton cm. scattering angle is 
9 = 90° (where Ui = ti), with eigenvalues that may be read off from Eq. 



In general, the eigenvalues of the anomalous dimension matrix, Eq. 



are 



Ai, 



Fn + F22±V(rii-r22)2 + 4Fi2F 
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and the eigenvectors are 



12 



Ai— Fii 
1 



(3.10) 



(3.11) 



for each eigenvalue Aj. These expressions will be useful when we discuss the 
diagonalization of the anomalous dimension matrices in Section 5. 



3.2 One- and two-loop expansions of the resummed 
cross section for qq -^ QQ 

We can expand the resummed heavy quark cross section to any fixed order 
in perturbation theory without having to diagonalize the soft anomalous 
dimension matrices. Such fixed-order expansions of resummed cross sections 
have also been done at one and two loops for the Drell-Yan process p^, Higgs 



production [Q , and heavy quark electroproduction [Q . We will present one- 
and two- loop expansions of the resummed cross sections for direct photon and 
W boson production in Section 10. 

First we give the one- loop expansion for qq -^ QQ. At one-loop the 
inverse Mellin transforms are trivial. The result is proportional to the Born 
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cross section, a^- 



M?) 



M) 



., ^ -^QQi ^^d the one-loop contributions from gf' and §2 , 
Eq. (|2.37|) , and Rer22, the real part of r22. In the DIS scheme the one-loop 
result is 



-.DIS(l) /. 2 + \ 



+ 



I- z 



-1 + 



+ Ca 



a.. 



a 



qq~^QQ' 



TT 



2C, 



Ml 



ln(l - z) 



C^ - + 8 In — - 2 - 2 ReLfl + 2 In — 



-"'^)-----(S' 




while in the MS scheme 



^MS(l) /, 2 + N 

(^qg^QQi^ - z,m ,S,tl,Ul) 



B 



""m^QQ^ 4Ci. 



ln(l-^) 

I- z 




C^(81n(— ") -2-2ReL^ + 21n(4 



Ui 



' im?s 



-3 In ( — + Re Lfl - In , 

ti / \tiUi 




. (3.13) 



Note that these expressions for the cross section are for fixed values of 
the QQ invariant mass. There are approximate one-loop results for heavy 
quark production in the literature for a single-particle inclusive cross section 
where the singular behavior is given in terms of the variable S4 = s + ti -|- mi 
pi]] . Therefore there are small differences in the results for the two cross 
sections due to the different phase spaces. Nevertheless, the cross sections 
are kinematically equivalent in the limit /3 — ^ 0, where S4 = 2m^{l — z). 

If one uses single-particle inclusive kinematics in the resummation proce- 
dure pi|, then one reproduces the results in Ref. |pT]] exactly. In addition one 



can add the Coulomb corrections to the one-loop expansions ||T9|, ^ (these 
corrections also appear in Ref. |^). Numerically the one-loop expansions of 
the NLL resummed cross section are very good approximations to the exact 



NLO cross sections at both the partonic and hadronic levels [Q, |3^ . 

One can expand the resummed cross section to higher orders and thus 
make predictions of perturbation theory at these higher orders near threshold. 
The inversion back to momentum space at any fixed order is straightforward. 
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For the two-loop expansion of the resummed cross section in the MS scheme 
and with /i = Q we have 



a 



MS(2)_ 
qq^QQ 



(1 - Z,m^,S,ti,Ui) 

\n^(l-z) 



R 



%-->QQ^ \ ^^F 



In^fl - z) 



+ 



1-z 



Cf 



-/5o + 12Cf (4 In (—] -ReLf^- l] 



Ml 



+ 6Ca -31n -i -In 



m?s 



+ ReL 



/3 



(3.14) 



Analogous results have been obtained in the DIS scheme, and also in single- 
particle inclusive kinematics in both schemes W^ - 



4 Soft anomalous dimension matrix for 
the process gg -^ QQ 

4.1 Eikonal diagrams and Ts for gg -^ QQ 

We continue with the soft anomalous dimension matrix for heavy quark pro- 
duction through gluon fusion, 



giPa, Ta) + g{pb, n) -> Q{pi, n) + Q{p2, r2) , 



(4.1) 



with momenta and colors labelled by the pj's and Tj's, respectively. Here r5 
is a 3 X 3 matrix. We use the same Mandelstam invariants as in the previous 
section, Eq. ( p.2|) . The UV divergent graphs are the same as in Fig. 2, where 
now the incoming eikonal lines represent gluons. In our calculations we use 
the eikonal rules as shown in Fig. 4 (see also the previous section). The 
gluon eikonal vertex is 

-gsr'%^A, (4.2) 

where we read the color indices a, b, c anticlockwise, and where A = +1(— 1) 
for the gluon located below (above) the eikonal line. 
We make the following choice for the color basis: 



Ci OraTi, 0r2rii 



C2 = d"^'"^ (T^; 



r'2r-u 



C, = ifranc^j.c^^ 



T2r\ 1 



(4.3) 
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-i v^/(-v • k+i 8) 



iv^/(-v • k+i 8) 



-v7(v-k+i8) 




v^/(vk+i8) 



Fig. 4. Feynman rules for eikonal lines representing gluons, in the process 
99 ~^ QQ- The gluon momentum flows out of the eikonal lines. Group 
matrices at the vertices on the left side are those of three-gluon vertices; 
those on the right are as for quark lines. 



where again the Tp are the generators of SU{3) in the fundamental repre- 
sentation, and d"'^'^ and f"''"^ are the totally symmetric and antisymmetric 
SU{3) invariant tensors, respectively. The anomalous dimension matrix in 
this color basis and in a general axial gauge is given by ||21|, |3T 



11 





131 
2 





r22 


— J- 31 


^31 


AN. ^ 31 


r22 



(4.4) 



where 
Til = 

r3i = 



r 



22 



TT 
TT 



-Cp{Lp + 1) + Ca[-- In (4i/,z/,) + 1 



TTZ 



In 



Ml 



TT L ^ 



In (4z/^z/fe) + 2 + In 



^1^1 
im?s 



+ Le 



m 



(4.5) 
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Here we haven't absorbed the function u^-'', Eq. ( |2.24| ), in the resuhs for 



Ts, as was done in Refs. pT| , |3T[| . Again we note that Ts is diagonahzed in 
this basis at absolute threshold, /5 = 0, and also for arbitrary /3, when the 
parton-parton cm. scattering angle is 6 = 90° , with eigenvalues that may 
be read off from Eq. (|4.5| ). 

The eigenvalues of Ts may be written as 



where 



Ai 

-^2,3 

X 



_(Xi/3_y + rn + 2r22), 



i(Xi/3 _Y)+ En + 2E22 ± ^tVsiX^/^ + Y) 



iru-r,,r-^{N'^-S)Tl,{T 



11 



22; 



+ 



3v^ 



256 



rL-(iVc-4)rL(rn-r22r 



+ o((iVc'-8)^-12(iV2-2))E^,(E 



11 



r22)^ 



nl/2 



and 



The eigenvectors of E5 are given by 



4) 



x-i/l 



2(A,-rii) 



N^r 



c^ 31 



4(Ai-r22) 

1 



for each eigenvalue Aj. 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



4.2 One- and two- loop expansions of the resummed 
cross section for gg -^ QQ 

Again, we may expand the resummed cross section for gg -^ QQ at one and 
two loops or higher. In this case the color decomposition is more compli- 
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cated I^T], ^]. The one-loop expansion in the MS scheme is 



a 



MS(1)_ 
gg^QQ 



[1- Z,m ,S,ti,Ui) 



a 



gg^QQ 



a. 



AC, 



TT 



ln(l - z) 



2C^ln 



fX^ 



1-2 



+ Ois^ggBqED 



^^^1 
+ T^^ 



1-z 

tiUi 



-Cf 



_ (t[+w| 



-CF + ^]ReLp 



m?s 



+ 



N^-l 



■m?s 



Nr. 






{CF~CA)ReLp 



1 In 



h 



(4.10) 



where a^ r,r, is the Born cross section, 

ti Ml 4m^s 

-DQED = — ^ :r + ~;^ — 

Ui ti tiUi 



1 - 



rn^s 



tlMi 



(4.11) 



and -ft", 



99 



{N^ 



is a color average factor. 



As we explained in the previous section, we cannot compare our one-loop 



expansion directly with the approximate NLO results of Ref. ||6T|, but as 
/3 — > our expression becomes identical to the /5 — *> limit of the results in 
[ |61[] . Again, we note that even for /? > our one- loop expansion is nearly 
the same as in Ref. [^]. If one uses single-particle inclusive kinematics in 
the resummation procedure then the agreement with Ref. [^ is exact |^ . 
The two-loop expansion of the resummed cross section in the MS scheme 
and with fi = Q is 



^MS(2) ,-, 2 J. ^ 



a 



0"„ 



99- 



' > 8Cl 



TT^ 



In^d 



ri„2 



-PoCa 



ln"(l 



7C 



ri„2 



ggJ^QED 



ln"(l - z) 

1-Z 



CaKnI 



\-z 

1-2 I „,2\ 



(tj + ul) 
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X 



N; In (^) - 2/V, (Cp - ^ ) Rei, - 2N^C, 



Analogous results may also be obtained in single-particle inclusive kinemat- 




ics 37 . 



5 Diagonalization procedure and 
numerical results 

As we saw in the previous two sections, the soft anomalous dimension ma- 
trices, Fg, are in general not diagonal. They are only diagonal at absolute 
threshold, /5 — >• 0, and at a scattering angle 9 = 90°. In these cases the 
exponentiated cross section has a simpler form. Numerical studies of the 
resummed cross section at ^ = 90° for top quark production at the Fermi- 
lab Tevatron and bottom quark production at HERA-B were presented in 



Ref. 33 



In general, however, we must find new color bases where F5 is diagonal 
so that the resummed cross section can take the simpler form of Eq. (|2.35|) . 
In this section we outline the required diagonalization procedure and apply 
it to heavy quark production, giving some numerical results for top quark 
production at the Tevatron. 

5.1 General diagonalization procedure 

The Born cross section can be written in an arbitrary color basis as a product 
of the hard components and the color tensors: 

a'' = Hjjc\cj. (5.1) 

In a diagonal basis, c'j = ckRkii the Born cross section can be rewritten as 

a^ = H,j{R-^)^j,{d)y^{R-^)Ki. (5.2) 

where R is made from the eigenvectors of the anomalous dimension matrix, 
and we have used c/ = c'f^R~^\. Then the resummed cross section is given by 

a-^ = Hu iR~-')Uc){cAR-')Ki e^- , (5.3) 
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where the exponent Ekl takes contributions at NLL from the eigenvalues 
\k and Al, and is given for /x = Q by Eq. (|2.36| ). 

Note that for the original color bases that we have chosen for all the 
subprocessess, we have ciCj = for / 7^ J; thus, our expressions simplify. 
Then the Born cross section can be written as 



a^ = Y.HH\ci\ 



(5.4) 



and the resummed cross section is given by 






(5.5) 



In the next two subsections we give more detailed results for the partonic 
processes involved in heavy quark production. 



5.2 Diagonalization and numerical results for the 
process qq -^ QQ 

The Born cross section for a process with a 2 x 2 anomalous dimension matrix, 
r^, such as qq -^ QQ, can be written in a general orthogonal color basis as 



a^ = Hulcil'^ + H22\c2\ 



(5.6) 



The eigenvalues and eigenvectors of Ts are given by Eqs. (|3.10| ) and (|3.11 



respectively. Then if C = (ci C2) is the original color basis, the diagonal 
color basis is 

C'^ic[ c'2) = CR, (5.7) 

where the columns of R are the eigenvectors of r5: 



R 



ei 62 



Ai — Fii A2— Fii 
1 1 



(5.^ 



The diagonalized anomalous dimension matrix is then given by 



pdiag ^ j^-lY^R 



Ai 
A2 



(5.9) 
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To write down the Born cross section in the diagonal basis we use the inverse 
of Eq. ( ^.71) , i.e. C = C'R~^, where the inverse of the matrix R is 



R' 



(Ai-rn)(A2-r 



11 



ri2(A2 — Ai^ 



4 



A2 — Til 

ri2 



Ai 



(5.10) 



We can thus express the old color basis tensors in terms of the new basis. 

In addition, we note that the Born cross section for the heavy quark 
production process qq -^ QQ is pure octet exchange so that there are further 
simplifications in the basis C = (cgingiet Coctct)- The Born cross section is 
then 



a 






= H22\C2\\ (5.11) 

with C2 = Coctet- Now, Coctet cau be written in terms of the diagonal basis as 

(Ai-rn)c;-(A2-rn)c'2 



(Ai - A2) 



(5.12) 



Then using the above equation we can rewrite the Born cross section in terms 
of Ic'ip, |c2p,CiC2*, where 



\r' P 

l'-l,2l 



^^c ■■- 12 



Ai,2 — Til 



+ 



N^-1 



and 



r' r'* 



A^2p2 
^^c '- 12 






(5.13) 



(5.14) 



(Ai-rn)(A2-rn)* ' 4 

To write the resummed cross section in momentum space a prescription 
must be chosen to invert the moment space exponentiated cross section. 



Here we use the simplest method of Ref . |jT8[ ; of course the diagonalization 
procedure is general and any prescription for the inversion can be used. 
The resummed partonic cross section is then given by |3^ 



<7is,m') 



2 

E 

«,i=i 



dcosO 



s~2ms^l'^ 



'^(0) 



dS4^fqq^ij{S4, 9) 



da 



qq,i3\"i 



S, S4, 



dsA 



(5.15) 
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where s^ = s + ti + Ui. The da^^^^^s^ S4, 9)/ds4 are components of the differ- 



ential of the Born cross section which is given by 



daf^(s,S4,e) 



-'nalKqqNcCp- 



sa 



4 ^■^■' ./(s - 54)^ -4sm2 

2 o„»^2\/-| I „„„2 /i\ I /i„„^2/'-| „„„2 , 



ds, ' ""^ •- ^ 4s4 

X (3(s- Sif -8sm^)(l + cos^6') + Asm\l -cos' 
The function / is given at NLL by the exponential 

fqq.ij = eyip[Egq,ij] = eXp [Egg + Egg{Xi, Xj)] , 

where in the DIS scheme 



E. 



DIS 



1 Ay 

77 



1^0 



VA^ d^ v2Cf 



1 



TT 



■(«.(0 

3Cr 



-i"2«)*-)i-if-(^)}' 



(5.16) 



(5.17) 



(5.18) 



with ujq = S4/(2m^) and A the QCD scale parameter. The color-dependent 
contribution in the exponent is 



Eqq{Xi, Xj 



"4U 

a;o CO I 



,/2/l2' 



a.o 



A2 



+ A' 



/2n2' 



«.« 



A2 



.^ 



(5.19) 

in both mass factorization schemes, where i = 1, 2. Here A' = A — ^^^^2 (see 
Eq. ( |2.37| )) where the A's are the eigenvalues of the soft anomalous dimension 
matrix. The cutoff Scut in Eq. ( ^.15| ) regulates the divergence of a^ at low S4: 

Scut > "S4^mm = 2?7i A/Q. 

After some algebra, the resummed partonic cross section at NLL accuracy 
is 
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crlfis,m') 



2 

E 



1 |.s-2msi/2 -j^ 

dcosO / ds4, -— — 

-1 "'Scut Al — A2 



d(T''^{s,S4,^ 



dsA 



X 



' ^^' r^ + lA r 



ip e^M'" + 



']^rL + |A,-rnr)e--- 



m: 



iV2 
c 



^r^^Re (e^-.-) - 2Re ((Ai - rn)(A2 - rn)*e^«^^- 



(5.20) 
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The explicit expressions for the quantities in the above equation are long 
but their derivation is straightforward. Numerical results for the partonic 
cross sections (resummed and one- loop expansions) are given in [Q. The 
one-loop expansion of the resummed cross section is a good approximation 
to the exact NLO result. 

The NLL resummed hadronic cross section is given by the convolution of 
parton distributions (pi/h, for parton i in hadron /i, with the partonic cross 
section 



(5.21) 



where a^^^TS^rn?) is defined in Eq. ( ^.20|) and tq = {m + \/m? -|- Scut)^/5'. 



Our numerical results for the tt production cross section at the Fermilab 
Tevatron with \fS = 1.8 TeV are shown in Fig. 5 as functions of the top 
quark mass. We use the CTEQ 4D DIS parton densities [^ |6^. Since 
the parton densities are only available at fixed order, the application to a 
resummed cross section introduces some uncertainty. The NLO exact cross 
sections, including the factorization scale dependence, are shown in Fig. 5 
along with the NLO approximate cross section, i.e. the one-loop expansion 
of the resummed cross section, calculated with Scut = and fi'^ = rri^. We 
note the excellent agreement between the NLO exact and approximate cross 
sections. 

As for the NLL resummed cross section, we note that its scale dependence 
is significantly reduced relative to that of the NLO cross section. In order to 
match our results to the exact NLO cross section we define the NLL improved 
cross section 



a, 



imp res NLO,approx |^ NLO, exact /r ooN 

''qq,ha.d ~ "^gg.had ~ "^qg.had + "^qg.had ; (^O.Z/J 

with the same cut applied to the NLO approximate and the NLL resummed 
cross sections. 

In Fig. 5 the hadronic improved cross section is shown for /i^ = m'^ along 
with the variation with Scnt- The variation of the improved cross section 
with change of cutoff is small over the range 30s4,min < Scut < 40s4,min- At 
m = 175 GeV/c^ and -vS = 1.8 TeV, the value of the improved cross section 
for qq —>■ ti with Scut/(2m^) = 0.04 is 6.0 pb compared to a NLO cross 
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Fig. 5. The NLO exact and approximate and the NLL improved hadronic ti 
production cross sections in the qq channel and the DIS scheme are given as 
functions of top quark mass for pp coUisions at the Tevatron energy, VS = 



1.8 TeV. The NLO exact cross section is given for n 



m 



(sohd curve), 



Am (lower-dashed) and m /A (upper-dashed). The NLO approximate cross 



section with s, 



cut 



is shown for /i 



m 



(lower dot-dashed). The NLL 



improved cross section, Eq. (|5.22| ), is given for Scut = 35s4,min (upper dot- 
dashed), 30s4^min (upper-dotted) and 40s4^min (lower-dotted). 
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section of 4.5 pb at /i = m. At the upgraded Tevatron with v>S' = 2 TeV, 
the corresponding value is 7.8 pb compared to a NLO cross section of 5.9 pb 
at /i = m. We find that the corrections relative to NLO are larger for larger 
scales. The gg channel is more complicated and will be discussed in the next 
subsection. Adding the gg contribution we predict a total cross section of 7 
pb at y/S = 1.8 TeV, in good agreement with experimental values from CDF, 

and DO, a^ = 5.5 ± 1.8 pb 



(J, 



7.Qtli pb 
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Gluon fusion is more 
important for 6-quark production at HERA-B where threshold resummation 
is also of importance [0, |2( 



5.3 Diagonalization for the process gg -^ QQ 

The Born cross section for a process with a 3 x 3 anomalous dimension matrix, 
r^, such as gg -^ QQ, can be written in a general orthogonal color basis as 



a 



R 



iflllCip + H22\C2? + -f^Ssksl 



(5.23) 



The eigenvalues and eigenvectors of Vs are given by Eqs. ( |4.(j| ) and (|4.9|), 

respectively. Then if C = (ci C2 C3) is the original color basis, the diagonal 
color basis is 

C'^{d, 4 d,)=CR, (5.24) 

where 



R 



ei 62 63] 



\e\ 


e\ 


e\] 




e\ 


el 


ei 
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1 





_aL. 



2(Ai-rii) 
4(Ai-r22) 

1 



Lar. 



2(A2-rii) 

4(A2-r22) 

1 



_ai. 



2(A3-rii) 

4(A3-r22) 

1 



(5.25) 
where, to simplify the expressions below, we denote by el the jth element 
of the zth eigenvector. The diagonalized anomalous dimension matrix is 

pdiag ^ ^-ip^^ 

To write down the Born cross section in the diagonal basis we use the 
inverse of Eq. ( |5.24| ), i.e. C = C'R~^, where the inverse of the matrix R is 



R- 



D 



-1 



r p2 _ 

^2 


-el 


e^ 


-e\ 


62^3 


-e\el 


2 
4- 


-el 


e\- 


-el 


2 1 
eie's- 


-e\e\ 
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-el 


e\- 


-e\ 


e;^e2 


-e\e\ 



(5.26) 
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with 

D = {e\- el)el - {e\ - e^)e^ - (e^ - e^)e? . (5.27) 

Then, the old basis may be expressed in terms of the new basis color tensors, 
and we can rewrite the Born cross section in terms of the diagonal basis as 



a 



B _ s:^ ^B 



E ^II, (5.2^ 



K,L=1 



r'* 



where a^^^ is the component of the Born cross section proportional to c'^^^^l 
The explicit expressions for these components are, again, long but they are 
straightforward to derive. Then the resummed cross section is given by 



^"- E ^IlS^-N (5.29) 



a 

K,L=1 



where the exponent Ej^l takes contributions at NLL from the eigenvalues 
\k and \l. 

6 Threshold resummation for dijet 
production 

In this section we review the resummation formalism relevant to the hadronic 
production of a pair of jets. We will follow the same methods as for heavy 
quark production but will encounter additional complications due to the 
presence of final state jets. 

6.1 Factorized dijet cross section 

We study dijet production in hadronic processes 

haiPa) + hipb) -^ Jlipu 6i) + J2(P2, ^2) + X (k) , (6.1) 

at fixed rapidity interval. 



1 In (P^ 

2 \Pi pI 



Ay = ^ln{^-^\, (6.2) 
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yjj = -.^A'-^^7^\- (6.3) 




with total rapidity, 



The jets are defined by cone angles 5i and 82. The introduction of cones 
removes all the final-state collinear singularities from the partonic cross sec- 
tion, which is then infrared safe, once the initial-state collinear singularities 
have been factored into universal parton distribution functions. We shall 
assume that the cones are small enough so that contributions proportional 
to 5j <C 1 may be neglected, but large enough so that as{Q) ln(l/5j) <^ 1, 
where Q is any of the hard scales of the cross section, such as the momentum 
transfer |^, |67[ . 



To construct the dijet cross sections, we define a large invariant, Mjj^ 
which is held fixed. A natural choice is the dijet invariant mass, 

M]j = {p,+P2)\ (6.4) 

which is the analog of Q^ for heavy quark production, but other choices are 
possible, for example the scalar product of the two jet momenta 

ML = 2pi-p2. (6.5) 

In both cases, large Mjj at fixed Ay implies a large momentum transfer 
in the partonic subprocess. The nature of the resummed cross section de- 
pends critically on this choice. As we shall see, the leading behavior of the 
resummed cross section is the same as for Drell-Yan or heavy quark produc- 
tion for the first choice for Mjj, Eq. (|6.4| ), while it is different for the second 
choice, Eq. (|6.5|). 

The dijet cross section is given in factorized form by 

dahaht,^.h.h{S, Mjj, y, Ay, 5i, 82) ^-^ f dxa dx 



dM]j dyjj dAy 



El wXfl CLXi, I / 2\ 

/ (Pfa/hAXa, 1^ ) 



fajb=q,g,9 ' 

.2\ Tj ( MJj IVljj 2n 



M^ M 

X (l^h/htixb, /i^) Hf^f^ I — ^, y, Ay, — ^, a,(/i^), 5i, 6; 



XaXi)0 fJ, 



2 



(6.6) 



where again H is the hard scattering and the 0's are parton distribution 
functions. The threshold for the partonic subprocess is given in terms of the 
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variable 2;, 

z = — ^ = —^ , (6.7) 

where, as before, S = {pa + VbY and s = XaXhS. At Zmax = 1 (partonic 
threshold) there is just enough partonic energy to produce the observed final 
state, with no additional radiation. The lower limit of z is 

z^i„ = r = -^ . (6.8) 

We may rewrite the cross section, introducing an explicit integration over 

z^ as 



dah^h.^j.j^jS, Mjj,y, Ay,6i, 62) _ y. f^ f dxg dxb 
dMjjdyjjdAy fa,h^q,q,g^'' ^ ^" ^^ 

^ (pfa/hA^a, f^ ) (j)f,/hSxb, ^^ ) S\z —I 6\^yjj--\n—j 

X Y. ^fah-*hf2{'^-z,— ^, Ay, asifi^), 61,62) , (6.9) 



fij2=q,q,g 



/i 



where again we have used the observation |^ that we can treat the total 
rapidity of the dijets as a constant, equal to its value at threshold. Thus we 
now have a simplified hard scattering function, o"/„/j-»/i/2- 

To calculate (yfah^hf2 at any order in perturbation theory, we construct 
the partonic cross section for the process /a + /b — > /i + /2, as above, 

daf^f^^j^j^{S,Mjj,Ay,6i,62) _ [^ f dxa dxb _^ ^^ 2^ 



dz / <pf^if^{Xa,ll ) 



dMjj dAy Jr J Xa Xb 

( M^ \ 

X (f>fa/fa{Xa, /X^) (l^h/hiXb, H^) 6\Z ^ j 

X Y. ^fah^hf2\^- z,^^, Ay, ols{^^'^), 61,62] , (6.10) 



fij2=q,q,g 



/^ 



where we have integrated over the total rapidity. We then factorize the 
initial-state coUinear divergences into the light-cone distribution functions 
0///, expanded to the same order in as, and thus obtain the perturbative 
expansion for the infrared-safe hard scattering function, a. As for heavy 
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Fig. 6. Refactorization for dijet production. The soft gluon function is as in 
Fig. 1(b). 



quark production, again we note that the leading power as 2; — ^ 1 conies 
entirely from flavor diagonal parton distributions. Moreover, we may sum, 
at leading power, over the flavors of the final-state partons /i, /2- 

By taking a Mellin transform of the rapidity-integrated partonic cross 
section (|6.10|) with respect to r, the above convolution becomes a product 



dr T 



N-l 



d'^fah^JiJ2iS, Mjj, Ay, Si, 62 



dMjj dAy 



x4>h/hiN,fi' 



f 

ai{N, Mjj/fx, a,(yu^), 5i, 82) , 



(6.11) 



where f denotes the partonic processes fa + fb—^fi + /2, and with d-{{N) = 
Jq dz z^~^(3"f(z), and (i){N) = J^ dx x^~^(j){x), as before. 

We can refactorize the cross section, as shown in Fig. 6, into hard compo- 
nents Hji, which describe the truly short- distance hard-scattering, center-of- 
mass distributions ip, associated with gluons collinear to the incoming par- 
tons, a soft gluon function Sli associated with soft gluons, and jet functions 
Jj, associated with gluons collinear to the outgoing jets. As before, I and 
L are color indices that describe the color structure of the hard scattering. 
The refactorized cross section may then be written as 



dr 



.N-l 



d(^fafb^JiJ2iS, Mjj, Ay, Si, 62) 



dMjj dAy 
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E E < -^^ Ay, a.(/.^) ) S^l> ( 9^, Ay, «,(/. 




+ 0{l/N). (6.12) 

This refactorization is similar to the heavy quark case except that now 
we have to include in addition outgoing jet functions in order to absorb the 
final state collinear singularities (the mass of the heavy quarks eliminates 
final state collinear singularities in heavy quark production). Expressions for 
the hard components Hji and the center-of-mass distributions ip were given 
in Section 2 in the context of heavy quark production. 

The soft function again represents the coupling of soft gluons to the par- 
tons in the hard scattering; this coupling, as we saw in Section 2.2, is de- 
scribed by eikonal lines. The construction of the eikonal cross section for 
dijet production is similar to the heavy quark case, apart from the outgoing 
jets. The collinear dynamics of the eikonal final-state jets are summarized 
by the matrix elements [^ 

where the ordered exponentials $^ are defined in Eq. ( |2.13|) , and with i = 1,2 



X (0| Trj f [<|.J{')t(oo,0;0)]|0(e|T[$J,f (oo,0;0)] j |0) , (6.13) 



and ^ a set of intermediate states which contribute to the weight Wi. 

We then construct moments of the soft function by dividing the moments 
of the eikonal cross section (defined in analogy to Eq. ( |2.15| ), see Ref. p8[] ) 
by the product of moments of the eikonal jets, Eqs. ( p.lTj ) and (|6.13|) , 



^/^ V jili^' {^, a.(/i^), e) M^^ {^, a.i^.^), 



X 



(6.14) 
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6.2 Resummed dijet cross section 



Comparing Eqs. ( |6.11| ) and ( |6.12| ), the refactorized expression for the MeUin 
transform of the hard scattering function is 



af{N) 






xT.Hni^,Ay,a.i,^) 



IL 
X J I 



/^ 



S'PJ^, Ay, a.(p^) 



(/l) 



A^, —^, a.(/i^), Si J(^f2) N, —^, as{fi'^),62 



/i 



/i 



(6.15) 



We discussed resummation for the ratio ip/(f) and the soft gluon function in 
the context of heavy quark production in Section 2. Thus, all we need to write 
down the resummed dijet cross section are expressions for the resummation 
of the final state jets. 



The moments of the final-state jet with Mjj = 2pi ■ p2 are given by |38 



Ji 



if.) 



'n.M" ^'-^^ 



i^ 



asin ),Si 



exp 



E(V^)(iV,M 



jj 



(6.16) 



with 



E[fAN,M^ 



JJ) 



,,N-1 



dz- 



{!--) dX 



1 — 2 [J(i-zy^ A 
+ B[f^[as{{l-z)Mlj)]} 



A(^) la^iXMh) 



(6.17) 



where the function A^^' is the same as in Eq. ( p. 21 ) and the lowest-order term 
in 5/ TN may be read off from the one- loop jet function. The results include a 
gauge dependence, which cancels against a corresponding dependence in the 
soft anomalous dimension matrix. The leading logarithms for final-state jets 
with Mjj = 2pi ■ p2 are negative and give a suppression to the cross section, 
in contrast to the initial-state leading-log contributions. 

The moments of the final-state jet with Mjj = {pi + ^2)^ are given by 
Eq. (§3|) with 



J fj, jJj 



(6.18) 
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where the first term in the series for C',rJas) may be read off from a one-loop 
calculation. The leading logarithmic behavior of the cross section in this case 
is not affected by the final state jets, so we always have an enhancement of 
the cross section at leading logarithm, as is the case for Drell-Yan and heavy 
quark cross sections. 

Using Eqs. ( |6.15|) , ( |2.27| ), ( p.31| ), and ( |6.16| ), we can write the resummed 
dijet cross section in moment space as 



af(iV) = Rff) expi ^ [E^f^\N,Mjj) 



-2 



Mj 



dfi' 



7/.(«.(/^"))-7/,:/,:(iV,«s(/^'^)) 



X exp ^: E[f^^{N,Mjjj 



X 



Tr H^'^ 



X Pexp 
X Pexp 



(— ^,A?/,a,(/i2)j 

"■'■''"^r«t(„.(,/^ 

^' -p(f) f i ■ 






12- 



(6.19) 



where E'^f^'^ is given by Eq. (|2.2CI| ), and E',f,-. is given by Eq. ( |6.17|) or Eq. (|6.18| ). 
This expression is similar to Eq. ( |2.33| ) for heavy quark production except 
for the addition of the exponents for the final-state jets. 

We give explicit expressions for the soft anomalous dimensions for the 
partonic processes relevant to dijet production in the next three sections. 



7 Soft anomalous dimension matrices for 
processes involving quarks 

In this section we present results for the soft anomalous dimension matrices 
for the processes qq — >• gg, qq —>■ qq, and qq — >• qq. The matrices for all these 
processes are 2x2. Hence the general diagonalization procedure for these 
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processes follows along the same lines as we discussed at the beginning of 
Section 5.2. 

Since the results are somewhat lengthy, we first introduce some notation 
which will simplify the expressions for Ts- 

7.1 Notation 

We consider partonic processes /„ {pa, ra) + fb {Pb, n) -^ /i (Pi, ri) + /2 (p2, ^^2), 
where the rj are color labels, and the Pi are momenta. To facilitate the 
presentation of the results for Ts we introduce the notation 

T = \n(—\+'ni, U = \n(—\+Tii, (7.1) 



where 



{Pa+Pbf , t = {.Pa-pif , U={pa-P2f, (7.2) 



are the usual Mandelstam invariants. Concerning the choice of physical chan- 
nel s, t, or u, for the definition of the color basis, we note that we can choose 
any channel. For the processes qq -^ qq, qq — *• qq and qq — *• qq, qg — *• qg and 
qg —>■ qg, as well as the process gg —>■ gg we will use t-channel bases, which 
seem to be the natural choice when analyzing forward scattering [Q. The 



processes qq -^ gg and gg -^ qq are better described in terms of s-channel 
color structures, so we will give results for them in s-channel color bases. 

Since the full cross section is gauge independent, i.e. independent of 
the choice of the axial gauge-fixing vector n'^, the gauge dependence in the 
product of the hard and soft functions, H^^^ S^^\ must cancel the gauge de- 
pendence of the incoming and outgoing jets, ip and J(/.). Now, the jets are 
incoherent relative to the hard and soft functions, so the gauge dependence of 
the anomalous dimension matrices Tg must be proportional to the identity 
matrix. Then we can rewrite the anomalous dimension matrix as 

{T^^)KL = {r^s')KL + SKL— E C(/,:)^(-lni/.-ln2+l-7r0, (7.3) 

with Cf^ = Cp {Ca) for a quark (gluon), and with Vi = (t>j-n)^/|np, as in Eq. 
(|3.9| ). Here the dimensionless and lightlike velocity vectors ff are defined by 

Pi = —7=-Vi , i = a,b,l,2, (7.4) 
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and satisfy vf = 0. 

For the subprocesses involved in dijet production we will present the 
explicit expressions for T^! . The full anomalous dimension matrix can then 
be retrieved from Eq. ( [7.3|) . 

We now give the results for the anomalous dimension matrices T^l for 
partonic processes involving quarks. 



7.2 Soft anomalous dimension for qq -^ qq 

First, we present the 2x2 soft anomalous dimension matrix for the process 



q (Pa, Ta) + q (Pb, Tb) ^ q (pi, Ti) + q (p2, ^2) 
in the t-channel singlet-octet color basis 

C2 = {''■ f) rirai-'- f) rtr2 ' 



(7.5) 



1 



1 



9N '''o-'^^ ''b''"i """ r)^rari,^rir2- 



(7.6) 



We find the soft anomalous dimension matrix 1^, |3S 



n 



2CfT -^U 



-2U 



JV, 



UT-2U) 



(7.7) 



The dependence on the logarithmic ratio T is diagonal in this t-channel color 
basis. In the forward region of the partonic scattering (T -^ —00), where F^ 
becomes diagonal, color singlet exchange is exponentially enhanced relatively 
to color octet |6^ . The general forms of the eigenvalues and eigenvectors of 
Ts are given by Eqs. ( p.lO|) and (|3.11|) ; explicit expressions are given in 
Ref. H. 

Finally we note that this result for F5 is consistent with the massless limit 
of the anomalous dimension matrix for the heavy quark production process 

gg nil 



qq 



7.3 Soft anomalous dimension for qq -^ qq and qq -^ qq 

Next, we consider the process 

q {Pa, Ta) + q (Pb, Tb) -^ q (pi, n) + q (pa, ^2) , (7.8) 
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in the t-channel singlet-octet color basis 

Cl — [-'- FJriraK-'- F)r2rb ~ ~ n ]\t ^rari^rf,r2 + Q"fa>"2"r(,ri ) 

C2 = 5rarAbr2- (7-9) 



The anomalous dimension matrix is given by [q^, |3^ 

-^^{T + U) + 2CfU 2U 
^U 2CfT 



J. as 
I s' — — 



(7.10) 



Again, the dependence on the logarithmic ratio T is diagonal in this t-channel 
color basis, and the color singlet dominates the octet in the forward region 
of the partonic scattering (T —^ — oo), where Vs becomes diagonal. 

Note that the same anomalous dimension matrix applies to the process 

Q{Pi,ri) + q{p2,r2) -^ q{Pa,ra) +q{Pb,rb) . (7.11) 

Again, the general forms of the eigenvalues and eigenvectors of Ts are 
given by Eqs. (|3.10|) and (|3.11| ). Explicit expressions for the eigenvalues and 



eigenvectors are given in Ref. |]39| . 

8 Soft anomalous dimension matrices for 
processes involving quarks and gluons 

Here we present the soft anomalous dimension matrices for the processes 
qq ~* 99^ 99 ^ qqi q9 ~^ q9i ^'^^ W ~~^ q9- The matrices for all these 
processes are 3x3. Moreover, we shall see that they are of the same general 
form as the anomalous dimension matrix for the heavy quark production 
process gg —>■ QQ. Hence the general diagonalization procedure for these 
processes follows the same lines as we discussed in Section 5.3 in the context 
of heavy quark production in the channel gg -^ QQ. 

8.1 Soft anomalous dimension for qq -^ gg and gg -^ qq 

First, we present the soft anomalous dimension matrix for the process 

q (Pa, ra) + q (Pb, n) -> g (pu ri) + g (pa, ^2) , (8.1) 
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in the s-channel color basis 



TbTa 



(8.2) 



We find m 



v., 



vr 







U -T 





^ (T + [/) ^f{U-T) 



2(t/-r) 



2Afc 



(8.3) 



The same anomalous dimension describes also the time-reversed process pl| , 

fi'(pi,ri) +fi'(p2,r2) ^ qiPa,ra) + q{pb,n). (8.4) 

We note that the anomalous dimension matrix, Eq. ( p.3|) , is of the general 
form of Eq. ([4.4|) . Then the general forms of its eigenvalues and eigenvectors 
are given by Eqs. ( ^.6| ) and ( [4.9| ); explicit expressions are given in Ref. p9[ . 
Finally we note that this result for Ts is consistent with the massless limit 
of the anomalous dimension matrix for the heavy quark production process 
gg^QQ [110. 



8.2 Soft anomalous dimension for qg -^ qg and qg 

Next, we consider the "Compton" process 

q (Pa, Ta) + g (Pb, Tb) "> q (pi, n) + g {p2, ra) . 
The soft anomalous dimension matrix in the t-channel color basis 

is given by |^ 

{Cp + Ca)T U 



qg 



(8.5) 



(8.6) 



r,., 



71 




2U 



CpT + ^U 



7V2-4] 



2Nc 



Cajj 
2 ^ 

CpT + ^U 



•7) 
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which also apphes to the process 

QiPuri) +g{p2,r2) ^ qiPa,ra) +g{.Ph,n) . 



(8.8) 



The dependence on the logarithmic ratio T is diagonal in this t-channel 
color basis, and the t-channel color singlet dominates in the forward region 
(T -^ — oo) where Ts becomes diagonal. Again, we note that the anomalous 
dimension matrix, Eq. ( ^.7|) , is of the general form of Eq. ( |4.4|) , so the general 
forms of its eigenvalues and eigenvectors are given by Eqs. ( ^.6| ) and ( [4.9| ). 
Explicit expressions for the eigenvalues and eigenvectors are given in Ref . [p^] . 



9 Soft anomalous dimension matrix and 
diagonalization for gg -^ gg 

9.1 Soft anomalous dimension matrix for gg -^ gg 

The final, and by far more complicated, process that we consider is 

g (Pa, ra) + g (Pb, n) -> gipi, ri) + g (pa, ra) . (9.1) 

The choice of a color basis, in which a four-gluon diagram can be expanded 



|69| , [rO|] , is more difficult in this case. In Ref. |3^ an initial overcomplete color 
basis of nine elements was found convenient for the calculations, resulting in 
a 9 X 9 anomalous dimension matrix. This was then reduced in a complete 
color basis to an 8 x 8 matrix. 

A complete color basis for the process gg -^ gg is given by the eight color 
structures |3^ 



Cl 
C5 



t 

4 
i r 

4 
i 

4 



frarbl^rir2l , ffartl irir2l 
frariljTf,r2l i ^faril £r,,r2l 



Pi{ra,rb]ri,r2) 



^(Jrari'^ri,r2 ? 



P8s{ra,n;n,r2) = -d''^'-''d'-'''' 



45 



C8 



PiiAi.'^a,H]ri,r2) 



_ fraTic rri,r2C 
O 



Pl0+wira,rb;ri,r2) = -((^r^r^'^r 



Or^r-i^rijri, 



_ rraric r 



ri,r2C 



1 1 

5 



(9.2) 



where the P's are t-channel projectors of irreducible representations of S'f/(3) 
[ [n]] , and we have used exphcitly A^c = 3. 



The soft anomalous dimension matrix in this basis is 1139 



TsxS 03x5 
05x3 r5x5 



with 



and 



r 



3x3 



TT 



3T 

3f/ 

3 (T + t/) 



(9.3) 



(9.4) 



r _ "^ 
J- 5x5 — 

TT 



6T 
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3r + f 
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3t/ 
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3r + f 





6t/ 
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-w 


3t/ 

4t/ 








EL 

4 
5 



^ _^ -2r + 4f/ 



(9.5) 



We note that the dependence on T is diagonal and in the forward region of 
the partonic scattering, T — > — oo, where T s becomes diagonal, color singlet 
exchange dominates. This has been a general trend for Ts for all the processes 
we analyzed in t-channel bases; suppression increases with the dimension of 
the exchanged color representation. 

The eigenvalues of the anomalous dimension matrix, Eq. (|9.3|), are 



Ai 

^7,8 



A4 = 

as 

2— 

TT 



3— T, 

TT 



At — As 



3^f/, 



Aa — Afi 



3— (T + f/), 

TT 



T + [/ T 2 VT2 -TU + W^ 



(9.6) 
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The eigenvectors have the general form 



0(5) 



, % i, Z, o , 6j 



0(3) 

J5) 



i = A...\ 



(9.7) 



where the superscripts refer to the dimension. The three-dimensional vectors 
ei are defined by 



,(3) 



5ii 

5i2 
5i3 



^=1,2,3, 



(9.8) 



while 0'-^'' and 0'-^-' are the five- and three-dimensional null vectors. The five- 
dimensional vectors e\\ 65 and Cq are given by p9[ 
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(9.9) 
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The expressions for e\ and eg are long but can be given succinctly by |^ 



,(5) 



" &i(AD " 


h2{K) 


h{K) 


b.iK) 


1 



(9.10) 



where 

and where the 6j's are given by 
3 



a.. 



(9.11) 



biiK) 



U^K 



-[80T^ + 103f/^ - 280UT^ - 300TU^ + AOAT'^U'^ 
+ (40T3 - 16f/^ - eOT^U + 52TU^)X 



2^^n 
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b2{K) 

HK) 



2K 



-[20T^ - 50UT + 44f/2 + (lOT - 5U)\[] 



-[AOT^ - 64f/3 - UOT'^U + ISOTt/^ 



+ (20T2 + 13f/^ - 20TU)\'i] , 

|;(2r+5t/-2A:), 



with 



(9.12) 

K' = 20T^ - 20UT + 21U^ . (9.13) 

With these resuhs for the eigenvalues and eigenvectors we have all the 
required elements for the diagonalization procedure. 



9.2 Diagonalization procedure for gg 
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We now give a discussion of the diagonalization procedure for the process 
99 ~^ 99- Since the anomalous dimension matrix is 8 x 8 its diagonalization is 
quite more complicated than for the other processes that we have considered. 
The Born cross section for gg -^ gg can be decomposed in the original 
basis, Eq. (|9.2|) , as 

a^ = Hn\ci\^ + H22\c2\^ + H^slc^l^ + H^^lcil^ 

+ i/sslcsr + Heelcel'' + Hjj\cj\^ + Hss\cs? . (9.14) 

Then if C = (ci C2 C3 C4 C5 cq c-j cs) is the original color basis, the diagonal 
color basis is 

C = (c^ C2 C3 C4 C5 Cg C7 Cg) = CR , 

where the columns of the matrix R are the eigenvectors of Y s''- 



(9.15) 



R = [ei 62 63 64 65 66 67 68] 



OsxS 



OsxS 
-^5x5 



(9.16) 



Here 



i?i 



5x5 



-15 

15 T 



-15 



6- i^- 

_15 T 
2 1/ 

3 
1 



2 




61 (A',) 61 (A- 



f + f^ 62(A'7) 62(A'g 



2 u 

15 T 
2 {/ 



3 T 
2U 



8J 

f-fij 63(A'7) 63(A'8) 
-3 64(A'7) 64(A'8) 
1 1 1 



(9.17) 
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-\diag 



and the diagonalized anomalous dimension matrix is T^f^ = R^^Ts'R- To 
decompose the Born cross section in the diagonal basis we use the relation 
C = C'R^^. The inverse of the matrix R is 



where 



R, 



5x5 





R-' = 


13X3 O3 
05x3 R, 


x5 
-1 
5x5 






12Ki 


iU^~bTU 
15Xi 


UT 
3X1 


t/2 

3i^i 


3t/2 





16t/2 
15^:2 





AU'^-STU 
3K2 


5X2 


C/2 


-U^+bTU 
15X3 


U'^-TU 
3X3 


C/2 
3X3 


31/2 


12X3 


20X3 


d.iX'r) 


d^m 


d,{\',) 


d,{K) 


rf5(A'7 


-diiX's) 


-d2{X's) 


-d,{\d 


-d. 


(A's) 


-d,{\ 



(9.1^ 



(9.19) 



f/^ 



f/-^ 



lOT^ + IST^f/ - 19Tf/2 + 24f/^ 
-A'(10T2-10Tf/ + 9f/2 



/(96iroir4) 



180T^ - 360T^U + 509T^U' - 339T^f/'' + 166Tf/* - 6f/ 



Here we have simplified the expression for the inverse of R by introducing 
the variables d, which are given as functions of the eigenvalues Aj, i = 7,8, 
by: 

diiXd 
d2{K) 

dsiK) 



- ^(60T^ 



U 



30T^ 
A' 



UOT^U + 163T2f/2 - 103Tf/^ 
- 4:5T^U + 59T^U^ - 32Tf/^ 4 



{- 48f/^) 
18f/^ 



/(12i^oi^5 



lOT^ - 15T^U + 13Tf/^ - 4f/^ 



/(12i^0^4) 



240T6 - eOOT^f/ + 972T^f/2 - SSST^'U^ + 583T2[/4 

- 205Tf/^ + 48f/6 - \i (aOT^ - lOOT^U + 152T^U^ 

- UST'^U^ + 70TU^ - 17U^)] /{2AK0K5) , 
2400T^ - 7200T^U + ISUOT^U^ - 15180T^f/^ + 12306T^U^ 

- 6321T^U^ + 2169TU^ - 264f/^ - A^ (400T^ - UOOT^U + 2U0T^U^ 

- 2280T^U^ + IQQQT^U^ - 72QTU^ + 181U^)] /{32K0K5) , (9.20) 



4(a: 
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where, again for brevity, we have introduced the notation 

Kq = VW^TUTTP, Ki = 4f/2 - AUT + ST^, 
K2 = AT^ - AUT + 5f/2, K3 = 5f/2 - QUT + ST^, 

'5 

-SOOT^f/ + lOOT^ (9.21) 



K^ = 20U^ - AATU^ + Q9T'^U^ - 50UT'^ + 25T\ 

K5 = lOOU^ - 300U^T + QOIU^T^ - 702U^T^ + QOIT'^U^ 



With these results we can now express the original basis in terms of 
the new diagonal basis color tensors. The Born cross section can then be 
rewritten in the diagonal basis as 



(7 



8 



E ^KL. (9-22) 



K,L=1 



where a^^^ is the component of the Born cross section proportional to dj^dl- 
While the calculation of these components is straightforward, the explicit 
expressions are quite long. 

The resummed cross section is then given by 



""'- E ^IlS^-N (9.23) 



a 

K,L=1 



where Ekl takes contributions at NLL from the eigenvalues \k and A^. 

10 Threshold resummation for direct photon 
and W boson production 

10.1 Factorization and resummed cross sections 

In this section we discuss resummation for direct photon and electroweak 
boson production. Threshold resummation may be important for these pro- 
cesses at large transverse momentum. Here we formulate the resummation 
in single-particle inclusive kinematics 1^1 . Resummation will follow, as we 
saw before, from the factorization properties of the cross section. 
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s 



(a) 




^^^^mr^ 




(b) 

Fig. 7. (a) Factorization for direct photon or W^ + jet production near 
partonic threshold, (b) The soft-gluon function S, in which the vertices c 
hnk ordered exponentials. 
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We consider hadronic processes of the form 

hA{pA) + hB{pB) ^ 7, W{Q) + X . (10.1) 

The partonic subprocesses contributing to direct photon oiW+ jet produc- 
tion are 

q{Pa) + g{pi,) — q{pj) + 7 , W{Q) (10.2) 

and 

q{Pa) + g(P6) ^ g{pj) + 7 , w^(Q) • (io.3) 

We define Mandelstam invariants 

S = {pa+Pb)\ t={pa-QY, U={pb-Q)\ (10.4) 

which satisfy s + t + u = Q'^ at threshold, with Q^ = M^, for W production 
and Q^ = for direct photon production. 

The factorized form of the cross section for direct photon or W + jet 
production is a convolution of the parton distribution functions with the 
hard scattering function, a: 

X a(s2,t,M,(5^a,(/i^)), (10.5) 

where we define S2 = s + t + u — Q'^, with Q = pw for W production and 
Q = P-y for direct photon production. The threshold region is given by S2 = 0. 
The cross section can be refactorized into distributions tp, defined in direct 
analogy with the cm. distributions that we presented in Section 2, hard 
components H = h*h, and a soft gluon function S. This refactorization, 
similar to the ones for heavy quark and dijet production, is shown in Fig. 7. 
There are some differences in the details of the formalism here as compared 
to the results in the previous sections because the resummation is now carried 



out in single-particle inclusive kinematics. More details are given in Ref. [^ 

The color structure of the hard scattering for direct photon or W boson 
production is much simpler than for heavy quark or dijet production. The 
color basis here consists of only one tensor; hence Ts is simply a 1 x 1 matrix, 
and no color traces or path ordering appear in the resummed expressions. 
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The resummation of the A^-dependence of each of the functions in the 
refactorized cross section leads to the expression in the MS factorization 
scheme in single-particle inclusive kinematics: 



d{N) 



expJ J2 [e^^'\N„p,-0 



-2 



2P«< dfl' 



7/.(a.(/i'^))-7/./.(iV.,«.(/^'')) 



X 






X S(l,(3,X,n,as{S/N^)) exp j J-2ReTs {as{^^!^] 

(10.6) 

where C,^ = p^/vS |^l| and /3j are the particle velocities. The first exponent 
E^^^^ has the same form as in Eq. ( |2.20 ) with A''^ = N{—u/s) and A*";, = 
N{—t/s). The exponent E',r x has the same form as in Eq. ( |6.17|) , with A^^^ 
given by ( p.21| ) and Birs given for quarks by |^ 



B' 



a. 



(9) 



c, 



TX 



and for gluons by 



B\ 



(9) 






^f 



16Ca 



-- + ln(2z/,) 



^ - 1 + ln(2z/,) 



(10.7) 



(10.8) 



10.2 Soft anomalous dimensions and expansions of 
the resummed cross sections 

Here we give the anomalous dimensions for direct photon and ly + jet pro- 
duction. The calculation follows the same lines as for heavy quark and dijet 
production. The main difference here is that the color basis consists of only 
one tensor, c = Tp, and there are three eikonal lines connecting at the color 
vertex. The one-loop eikonal vertex corrections for the partonic subprocesses 
in direct photon and W + jet production are shown in Fig. 8. 
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Fig. 8. One-loop eikonal vertex corrections for partonic subprocesses in 
direct photon and W + jet production. 



The anomalous dimension for the process q{pa) + g{Pb) -^ lipj) + liQ) 
or the process q{pa) + g{pb) -^ q{pj) + W{Q) is given by 



277 I 

+ C 



-u 



21n(^— j-ln(H„z/,J + 2 
In ( -") - ln(2z/g) + 1 - vri 



:i0.9) 



The one-loop expansion of the resummed cross section for direct-photon 
production is given in Ref. |^. The authors have found agreement with the 
exact NLO results in Ref. 0. Numerical results are given in Ref. [j72|. We 
may also expand the resummed cross section for W + jet production. The MS 
one- loop expansion at NLL of SijQ'^ in single-particle inclusive kinematics is 



-MS(1) 



%J,M/(s2,S,t,M,(5 



- {Cp + Ca) In f ^ 



(J, 
1 

S2 



B 

qg^qW 



a 



^{{Cf + 2Ca) 



TX 



Hs2lQ^ 



S2 



+ 



>^'H'€ 



1 

S2 



;io.io) 



where cr^^g^y denotes the Born cross section for this partonic channel. This 
result agrees with the exact NLO expressions in Ref. [0. 
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The anomalous dimension for the process q{pa) + <i{j>b) -^ giP-j) + liQ) 
or the process q{j)a) + <i{pb) — ^ gipj) + ^iQ) is given by 



2^ 



{Cf [- ln(4z/gZ/g-) + 2 - 27ri] 



+ C^ 



In ("^ ) - ln(2z/g) + 1 



TTi 



(10.11) 



Again, the one-loop expansion for direct-photon production is given in 
Ref. ^1 and agreement has been found with the exact NLO results in [^]. 
The MS one-loop expansion at NLL of S2/Q'^ for W + jet production in this 
channel and in single-particle inclusive kinematics is 



-MS{1) 



-MS 111 / , ,o2\ 



2C...^ 



+ Ca 



n 



f 



B 

qq^gW ' 

1 
11 



as 



TC 



{4Cp - Ca) 



In(s2/Q2) 



S2 



6Ca 12 



In 



2Cpln 



tu 



tu 



1 

S2 



(10.12) 



which again agrees with the exact NLO results in Ref. |0|. 

One may of course expand the resummed cross sections to two loops or 
higher for both direct photon []72| and W boson |^ production. For example, 
the MS two-loop expansion at NLL of 1 — tt; = S2/{s + 1) for direct photon 
production in the channel qq — > (77 is 



^MS(2), 



w, s,v) 



a 



R 



a: 



gq-^gi ^2 



'sCl - ACfCa + ^ 



In'fl 



w] 



w 



+ 



-12Cl In 



1-v 



..w^ +c, 



+ CfCa {9\n{l 
3 



v) — 31nt> + 



Ca 



-- n(l -v) ^H 

2 ^ ^ 4Ca 8 



— PoiCp ~ -Ca 
where v = 1 + t/s. 



rin^ri 



w] 



1 — w 



(10.13) 
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11 Conclusion 

We have reviewed the resummation of threshold logarithms for heavy quark, 
dijet, direct photon, and W boson production in hadronic collisions. Resum- 
mation follows from the factorization properties of the cross sections and is 
influenced by the color exchange in the hard scattering. We have presented 
full results for the soft anomalous dimension matrices for all the relevant par- 
tonic subprocesses. The one-loop expansions of the resummed cross sections 
agree with exact NLO calculations and one can expand the cross sections 
to two-loops or higher orders. We have discussed the general diagonaliza- 
tion procedure that can be implemented in the calculation of the resummed 
cross sections. Numerical results have been presented for top quark produc- 
tion at the Fermilab Tevatron. The resummation formalism is quite general 
and can be applied as well to the calculation of transverse momentum or 
other differential distributions for a variety of processes, such as heavy quark 
production. 
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A Evaluation of one-loop eikonal vertex 
corrections 

Here we give some calculational details for the one-loop corrections of Figs. 

2 and 8. 

In our calculations we use Feynman eikonal rules as discussed in Sections 

3 and 4, and a general axial gauge gluon propagator, 

D^'^ik) = T^N^'ik), N^^ik) = r-^^—T^+n'^^,. (A.l) 
k^ + le n ■ k [n ■ ky 
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with n^ the axial gauge-fixing vector. 

We denote the kinematic part of the one-loop vertex correction to c/, with 
the virtual gluon linking lines Vi and f j, as u!ij{6iVi, SjVj, Aj, Aj). The 6's and 
A's are defined as in Sections 3 and 4 except that here we use A = — i (+i) 
for a gluon located below (above) the eikonal line in order to present the 
results for both quarks and gluons in a uniform fashion. 

The expression for ujij is then 



uJij{SiVi,6jVj,Ai,Aj) = g"^ 
AiVi-n P 



Ai Aj Vi-Vj 



d^q —i j 
(27r)" q'^ + ie\ {SiVi-q + ie){6jVj-q + ie) 



AjVyn 



P 



{6iVi-q + ie) (n-q) {SjVj-q + ie) (n-q) 



+ n^ 



P 



[n-qY^ 



(A.2) 



with (yf^ = Alias, and where P stands for principal value, 



1 



(q-n)'^ 2 \{q-n + ie)(^ 
We may rewrite ( |A.2|) as 

uJij{6iVi,5jVj,Ai,Aj) = S, 



+ (-iy 



[—q ■ n + ieY 



(A.3) 



Ii{5iVi,5jVj) - -hikvi^n) - -l2{SiVi, -n) 



- -hiSjVj, n) - -hiSjVj, -n) + h{n^ 



(A.4) 



where Sij is an overall sign 

Sij = AiA,5i5j. (A.5) 

We now evaluate the ultraviolet poles of the integrals. For the integrals 



when both Vi and Vj refer to massive quarks we have (with e = 4 — ra) [^, ^ 



UV pole 
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rUV pole 
Jo 



rUV pole 
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(A.6) 
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where the Lp is the velocity-dependent eikonal function 

1 



'P 



2mVs A 1 - /5 ^ 

m :: + Til 



/3 



1+/9 



(A.7) 



with [3 = Jl — Am? / s. The Li and Lj are rather comphcated functions of 
the gauge vector n. Their contributions are cancelled by the inclusion of self 
energies. Their explicit expressions are: 



where 

LA±n) 



1 



Li = - [Li{+n) + Li{-n)] 



Wi ■ n\ 



(A.8) 
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+ In . 

5{-^n) r? — \vi ■ n\ + J {vi ■ riY — 2m?r? / s ^ 

with 5{n) = \vi ■ n\/{vi ■ n). Then, using Eq. ( |A.4[ ), we get 



a. 



Sij — [Lp + Lj + Lj — 1] . 
Tre 



(A.9) 



(A.IO) 



The contribution of the heavy quark self energy graphs in Fig. 2(b) to the 
anomalous dimension matrices for heavy quark production in either partonic 
channel is {as/Ti)CF{Li + L2 — 2)6ij, which cancels the gauge-dependent 
contribution from the Uu in Eq. ( |A.10|) . 

When Vi refers to a massive quark and Vj to a massless quark we have 
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where 



{va ■ nf 



\n\ 



(A.12) 



and Vij = Vi ■ Vj. Note that the double poles cancel in the sum over the J's 
and we get 



ujij{6iVi,5jVj,Ai,Aj 



vre 



2 \2m? 



^l^' 



+ Li + - In z/j - 1 



. (A.13) 



Finally, when both vi and Vj refer to massless quarks we have [p5| , |21| , pi 
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(A. 14) 



Again, note that the double poles cancel in the sum over the /'s and we get 

(A.15) 



ujij{6iVi, 6jVj, Aj, Aj) = Sij — 



a. 
ne 



In ( 6i 6j y) + 2 Hi^'^i^i) - 1 



In order to obtain contributions to the different entries of the matrix of 
renormalization constants, the above expression has still to be multiplied by 
the color decomposition of its corresponding diagram into the basis color 
structures ETI, IHl 
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